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Abstract. We present the new technique for the determination of the low-mass slope of the 
stellar mass function from the pixel-space fitting of integrated light spectra. This technique 
is the extension of the NBursts full spectral fitting technique. It can be used to constrain 
the stellar initial mass function (IMF) of compact stellar systems having high relaxation 
timescales (Hubble time or more). We used Monte-Carlo simulations with mock spectra to 
test the technique and conclude that: (1) age, metallicity, and IMF can be very precisely 
^H. determined in the first "unrestricted" variant of the code for high signal-to-noise datasets 

q (S/N= 1 00, R=7000 give the uncertainty of alpha of about 0. 1 ) ; (2) adding the M/L informa- 

5^ ' tion significantly improves the precision and reduces the degeneracies, however systematic 

, errors in M/L will translate into offsets in the IMF slope. 
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ly-} ■ Globular clusters (GC) and ultracompact approach for determination of parameterized 

' dwarf galaxies (UCD) do not have a dark mat- the line of-sight velocity distribution and 

■ ter according to the modern studies. Hence, star formation history. NBursts uses the pre- 
their dynamical masses derived from velocity computed grid of the high-resolution pegase.hr 
dispersion measurements equal to their stellar (Le Borgne 2004) simple stellar population 

(f) ■ masses computed from the stellar population (SSP) models. 

parameters and the slope of the stellar current in t he first version of our new approach 

mass function. If we use the spectra of compact an observed spectrum is fitted against an opti- 

stellar systems having high dynamical relax- ma i template represented by a linear combina- 

ation timescales (Hubble time or longer) then t i on G f sSPs each of them characterised by the 

the effects of the mass segregation and tidal age , metallicity and the low-mass IMF slope a, 

evaporation of low-mass stars can be neglected determined in the same minimization loop, us- 

and the current stellar mass function should be i ng 3-dimensional cubic spline interpolation on 

equivalent to IMF. the grid of high-resolution pegase.hr SSP mod- 

Our technique is the extension of the els. We deal with observed spectra of compact 

NBursts (Chilingarian 2007). NBursts is the stellar systems which can be well reprensented 
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Fig. 1. Monte Carlo simulations using mock spec- 
tra (R=7000) demonstrating the precision of the 
IMF slope determined by our technique as well as 
degeneracies between parameters and systematic er- 
rors cause by wrong estimates of (M/L) dyn . Red con- 
tours correspond to the results of the fitting by the 
first version of our technique. The blue contours 
correspond to the second version with the exter- 
nal imposed M/L value. Green and black contours 
show how the fitting results change when we use 
the external M/L value offset from the input one 
by -15% and +15% correspondingly. The degen- 
eracies between metallicity and age (not presented 
at this figure) are weak (Ar=0.5 Gyr approximately 
corresponds A[Fe/H]=0.05 dex). 



by single component SSP models. The model 
grid has 25 nodes in age (10 Myr to 20 Gyr), 
10 nodes in metallicity ([Fe/H] from -2.5 to 
1.0 dex) and 10 nodes in the low-mass IMF 
slope (0.4 to 3. 1 with a step of 0.3). 

The second version of the our tech- 
nique is mathematically identical to the orig- 
inal nbursts technique but the specific grid 
of input SSP models is constructed and sup- 
plied for every observed spectrum. Here we 
assume that the object being studied contains 
no dark matter. In this case, (M / L)& yn =(M / L)*. 
(M/L)d yn were determined from the analysis of 
their internal structure and observed velocity 
dispersion profiles and are available in the liter- 
ature (Mieske 2008; Chilingarian 2008, 2010). 
(M/L)„ can be derived from pegase.2 stellar 
population models for every set of (t, [Fe/H], 
a) and for every given t and [Fe/H] in the 



old stellar population regime this function is 
monotonic. Hence, if we know (M/L)d yn and 
impose the zero dark matter constraint, for ev- 
ery point on the (f, [Fe/H]) plane we can find 
the value of a such as (M/L)» in that point 
of the parameter space equals to (M/L)d yn - 
Thus, we can compute a grid of SSP models 
in the age-metallicity space varying the low- 
mass IMF slope so that the (M/L), values are 
constrant all over the grid and equal to the 
"external" (M/L)d yn - Along with this grid, we 
will also map the behaviour of a as a func- 
tion of t and [Fe/H]. Then, if we feed this SSP 
grid to the standard nbursts full spectral fit- 
ting procedure and determine the pair of best- 
fitting values of age and metallicity for some 
fixed (M/L), we will automatically get the a 
value corresponding to this best-fitting solution 
which will measure the low-mass IMF slope of 
the stellar population in the stellar system be- 
ing studied. 

We also applied our technique to observed 
intermediate and high-resolution GC and UCD 
spectra obtained with Gemini and VLT. High 
signal-to-noise spectra having R = 2000 at 
3900 < A < 6800 A (Francis 2012) yield to 
cr(a) = 0.2 . . . 0.4 and 0.1... 0.2 respectively 
in the first and second verions of our approach. 

We conclude that by applying our tech- 
nique to high-quality optical observations of 
extragalactic GCs and UCDs we are able to 
reach better precision of the IMF determina- 
tion than that made with direct star counts in 
nearby open clusters and check the IMF uni- 
versality hypothesis (Kroupa 2002). 
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